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PACS 96 . 60 . qd - Sun spots, solar cycles 

Abstract. - Solar activity is studied using a cluster analysis of the time- fluctuations of sunspot 
number. In the historic period (1850-1932) the cluster exponent a ~ 0.37 (strong clustering) for 
the high activity components of the solar cycles. In the modern period (last seven solar cycles: 
1933-2007) the cluster exponent was a ~ 0.50 (random, white noise-like). Comparing these 
results with the corresponding data from laboratory experiments on convection it is shown, that 
in the historic period emergence of sunspots in the solar photosphere was dominated by turbulent 
photospheric convection. In the modern period, this domination was broken by a new more active 
dynamics of the inner layers of the convection zone. Cluster properties of the solar wind magnetic 
field and the aa-geomagnetic-index also support this result. Long-range chaotic dynamics in the 
solar activity is briefly discussed. 



Introduction. — The sunspot number is the main di- 
rect and reliable source of information about the solar dy- 
namics for a historic period. In a recent paper [1], for 
instance, results of a reconstruction (based on radiocar- 
bon concentrations, sec also [2]) of the sunspot number 
were presented for the past 11,400 years. Analysis of this 
reconstruction showed an exceptional level of solar activity 
during the past seven solar cycles (1933-2007). Sunspots 
appear as the visible counterparts of magnetic flux tubes 
in the convective zone of the sun. Since a strong magnetic 
field is considered as a primary phenomenon that controls 
generation of the sunspots the crucial question is: Where 
has the magnetic field itself been generated? The location 
of the solar dynamo is the subject of vigorous discussions 
in recent years. A general consensus had been developed 
to consider the shear layer at the bottom of the convec- 
tion zone as the main source of the solar magnetic field [3] 
(see, for a recent review [4]). In recent years, however, the 
existence of a prominent radial shear layer near the top 
of the convection zone has become rather obvious and the 
problem again became actual. The presence of large-scale 
meandering flow fields (like jet streams), banded zonal 
flows and evolving meridional circulations together with 
intensive multiscale turbulence shows that the near sur- 
face layer is a very complex system, which can significantly 
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Fig. 1: Daily sunspot number (SSN) vs time [12]. The dashed 
straight line separates between historic and modern periods. 



affect the processes of the magnetic field and the sunspots 
generation. There could be two sources for the poloidal 
magnetic field: one near the bottom of the convection zone 
(or just below it [3]), another one resulting from an active- 
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Fig. 2: The standard deviation for 5n T vs r for the level 
SSN=85 (historic period) in log-log scales. The straight line 
(the best fit) indicates the scaling law Eq. (1). 



Fig. 3: The standard deviation for Sn T vs r for the level 
SSN=125 (modern period) in log-log scales. The straight line 
(the best fit) indicates the scaling law Eq. (1). 



region tilt near the surface of the convection zone. For the 
recently renewed Babcock-Lcighton [5], [6] solar dynamo 
scenario, for instance, a combination of the sources was 
assumed for predicting future solar activity levels [7], [8]. 
In this scenario the surface generated poloidal magnetic 
field is carried to the bottom of the convection zone by 
turbulent diffusion or by the meridional circulation. The 
toroidal magnetic field is produced from this poloidal field 
by differential rotation in the bottom shear layer. Destabi- 
lization and emergence of the toroidal fields (in the form 
of curved tubes) due to magnetic buoyancy can be con- 
sidered as a source of the pairs of sunspots of opposite 
polarity. The turbulent convection in the convection zone 
and. especially, in the near-surface layer captures the mag- 
netic flux tubes and cither disperses or pulls them trough 
the surface to become sunspots. 

The magnetic field plays a passive role in the photo- 
sphere and does not participate significantly in the turbu- 
lent photospheric energy transfer. On the other hand, the 
very complex and turbulent near-surface layer (including 
photosphere) can significantly affect the process of emer- 
gence of sunspot. The similarity of light elements prop- 
erties in the spot umbra and granulation is one of the 
indication of such phenomena. 

The paper reports a direct relation between the fluctua- 
tions in sunspot number and the temperature of turbulent 
convection in the photosphere for the historic period. This 
relation allows for certain conclusions about the genera- 
tion mechanisms of the magnetic fields and the sunspots. 
In the modern period (last seven solar cycles: 1933-2007, 
cf. [1], [2]) the relative role of the surface layer (photo- 
sphere) in the process of emergence of sunspots decreased 
in comparison with the historic period (the seven solar 
cycles preceding 1933), implying a drastic increase of the 
relative role of the inner layers of the convection zone in 
the modern period. 



Time-clustering of fluctuations. — In order to ex- 
tract a new information from sunspot number data we 
apply the fluctuation clustering analysis suggested in the 
Ref. [9]. For a time depending signal we count the num- 
ber of 'zero'-crossings of the signal (the points on the time 
axis where the signal is equal to zero) in a time interval r 
and consider their running density n T . Let us denote fluc- 
tuations of the running density as Sn T = n T — (n T ), where 
the brackets mean the average over long times. We are 
interested in scaling variation of the standard deviation of 
the running density fluctuations (Sri 2 ) 1 / 2 with r 
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For white noise signal it can be derived analytically [10], 
[11] that a = 1/2 (see also [9]). The same consideration 
can be applied not only to the 'zero'-crossing points but 
also to any level-crossing points of the signal. 

In Fig. 1 we show the daily sunspot number (SSN) for 
the period 1850-2007 [12]. Even by eye one can see that 
the modern period (last seven cycles, 1933-2007) is dif- 
ferent from the corresponding historic period: 1850-1933. 
Therefore, in order to calculate the cluster exponent (if 
exists) for this signal one should make this calculation 
separately for the modern and for the historic periods. 
We are interested in the active parts of the solar cycles. 
Therefore, for the historic period let us start from the 
level SSN=85. The set of the level-crossing points has a 
few large voids corresponding to the weak activity peri- 
ods. To make the set statistically stationary we will cut 
off these voids. The remaining set (about 10 4 data-points) 
exhibits good statistical stationarity that allows to calcu- 
late scaling exponents corresponding to this set. Fig. 2 
shows (in the log-log scales) dependence of the standard 
deviation of the running density fluctuations {Sn 2 ) 1 ^ 2 on 
r. The straight line is drawn in this figure to indicate the 
scaling (1). The slope of this straight line provides us with 
the cluster-exponent a = 0.37±0.02. This value turned to 
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Fig. 4: The standard deviation for Sn T vs r for the temperature 
fluctuations in the Rayleigh-Bernard convection laboratory ex- 
periment with Ra — 1.4 x 10 [14]. The straight line indicates 
the scaling law Eq. (1). 
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Fig. 5: The standard deviation for 5n T vs. r for radial compo- 
nent B r of the interplanetary magnetic field, as measured by 
the ACE magnetometers for the last solar cycle (hourly aver- 
age [17]). 



be insensitive to a reasonable variation of the SSN level. 
Results of analogous calculations performed for the mod- 
ern period are shown in Fig. 3 for the SSN level SSN=125. 
The calculations performed for the modern period provide 
us with the cluster-exponent a = 0.5±0.02 (and again this 
value turned to be insensitive to a reasonable variation of 
the SSN level). 

The exponent a ~ 0.5 (for the modern period) indicates 
a random (white noise like) situation. While the exponent 
a ~ 0.37 (for the historic period) indicates strong clus- 
tering. The question is: Where is this strong clustering 
coming from? It is shown in the paper [9] that signal pro- 
duced by turbulence exhibit strong clustering. Moreover, 
the cluster exponents for these signals depend on the tur- 
bulence intensity and they arc nonsensitive to the types 
of the boundary conditions. Fortunately, we have direct 
estimates of the value of the main parameter characteriz- 
ing intensity of the turbulent convection in photosphere: 
Rayleigh number Ra ~ 10 11 (see, for instance [13]). In 
Fig. 4 we show calculation of the cluster exponent for the 
temperature fluctuations in the classic Rayleigh-Bernard 
convection laboratory experiment for Ra ~ 10 11 (for a 
description of the experiment details see [14]). The cal- 
culated value of the cluster exponent a = 0.37 ± 0.01 co- 
incides with the value of the cluster-exponent obtained 
above for the sunspot number fluctuations for the historic 
period. If the value of the Rayleigh number Ra in the 
photosphere for the historic period has the same order 
as for the modern period: Ra ~ 10 11 (see next Section), 
then one can suggest that the clustering of the sunspot 
number fluctuations in the historic period is due to strong 
modulation of these fluctuations by the turbulent fluctu- 
ations of the temperature in the photospheric convection. 
This seems to be natural for the case when the photo- 
spheric convection determines the sunspot emergence in 
the photosphere. However, in the case when the effect 



of the photospheric convection on the SSN fluctuations is 
comparable with the effects of the inner convection zone 
layers on the SSN fluctuations the clustering should be 
randomized by the mixing of the sources, and the clus- 
ter exponent a ~ 0.5 (similar to the white noise signal). 
The last case apparently takes place for the modern pe- 
riod. Since the Rayleigh number Ra of the photospheric 
convection preserves its order Ra ~ 10 11 with transition 
from the historic period to the modern one (see next Sec- 
tion), we can assume that just significant changes of the 
dynamics of the inner layers of the convection zone (most 
probably - of the bottom layer) were the main reasons for 
the transition from the historic to the modern period. 

Magnetic fields in the solar wind and on the 
Earth. — Although in the modern period the turbulent 
convection in the photosphere have no decisive impact on 
the sunspots emergence, the large-scale properties of the 
magnetic field coming through the sunspots into the pho- 
tosphere and then to the interplanetary space (so-called 
solar wind) can be strongly affected by the photospheric 
motion. In order to be detected the characteristic scal- 
ing scales of this impact should be larger than the scaling 
scales of the interplanetary turbulence (cf. [15], [16]). In 
particular, one can expect that the cluster-exponent of the 
large-scale interplanetary magnetic field (if exists) should 
be close to a ~ 0.37. In Fig. 5 we show cluster-exponent 
of the large-scale fluctuations of the radial component 
B r of the interplanetary magnetic field. For computing 
this exponent, we have used the hourly averaged data ob- 
tained from Advanced Composition Explorer (ACE) satel- 
lite magnetometers for the last solar cycle [17]. As it was 
expected the cluster-exponent a ~ 0.37 ±0.02. Analogous 
result was obtained for other components of the interplan- 
etary magnetic field as well. 

At the Earth itself the solar wind induced activity is 
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Fig. 6: Cumulative aa-index vs. time. Daily aa-index was 
taken from [19]. The arrow indicates beginning of the transi- 
tional solar cycle. 



Fig. 7: The standard deviation for Sn T vs. r (in log-log scales) 
for the daily aa-index [19]: circles correspond to the historic 
period and crosses correspond to the modern period. 



measured by geomagnetic indexes such index (in 

units of 1 nT). This, the most widely used long-term ge- 
omagnetic index [18], presents long-term geomagnetic ac- 
tivity and it is produced using two observatories at nearly 
antipodal positions on the Earth's surface. The index is 
computed from the weighted average of the amplitude of 
the field variations at the two sites. We are interested in 
the clustering properties of the aa-index and their relation 
to the modulation produced by the photospheric convec- 
tive motion. The point is that the data for the aa-index 
are available for both the modern and the historic periods. 
This allows us to check the suggestion that the Rayleigh 
number of the photospheric convection has the same order 
for the both mentioned periods. The transition between 
the two periods is clear seen in Figure 6, where we show 
the cumulated aa-index a(t) 

Saa(t) = f a(t')dt' (2) 

Jo 

The arrow in this figure indicates beginning of the transi- 
tional solar cycle. 

Figure 7 shows cluster-exponents for both the modern 
and the historic periods calculated for the low intensity 
levels of the aa-index (for the historic period the level 
used in the calculations is aa-index=10nT, whereas for the 
modern period the level is aa-index=25nT, the aa-index 
was taken from [19]). The low levels of intensity were 
chosen in order to avoid effect of the extreme phenomena 
(magnetic storms and etc.). The straight line in Fig. 7 is 
drawn to show the expected value of the cluster-exponent 
a ~ 0.37 for the both periods. It means that indeed for 
both the historic and the modern periods the Rayleigh 
number Ra ~ 10 11 (cf. Fig. 4) in the solar photosphere. 

Chaotic sun. — The long-range reconstructions of the 
sunspot number fluctuations (see, for instance, Refs. [1], 



[2]) allow us to look on the solar transitional dynamics 
from a more general point of view. In figure 8 we show a 
spectrum of such reconstruction for the last 11,000 years 
(the data, used for computation of the spectrum, is avail- 
able at [20]). The spectrum was computed using the max- 
imum entropy method as in Rcf. [21] (see below). A semi- 
logarithmical representation was used in the figure to show 
an exponential law 

E(f) ~ e-ff" (3) 

The straight line is drawn in Fig. 8 to indicate the expo- 
nential law Eq. (3). Slope of the straight line provides us 
with the characteristic time scale r = 1/ f — 176y(±7y). 
The exponential decay of the spectrum excludes the pos- 
sibility of random behavior and indicates the chaotic be- 
havior of the time series. It is well known that low-order 
dynamic (deterministic) systems have as a rule exponen- 
tial decay of E(f) in the chaotic regimes (see, for in- 
stance, [21]). As for the delay-differential equations with 
chaotic attractors it is interesting to compare Fig. 8 with 
figure 3 of the Ref. [22]. It should be noted that the 176y 
period is the third doubling of the period 22y. The 22y 
period corresponds to the Sun's magnetic poles polarity 
switching. 

The exponential spectrum can be also produced by a 
series of Lorentzian pulses with the average width of the 
individual pulses equals to t (though, the distribution of 
widths of the pulses should be fairly narrow to result in 
the exponential spectrum). 

In Fig. 8 a local maximum corresponding to the fre- 
quency /o and its first harmonics have been indicated by 
arrows. It should be noted that the harmonic 2/ corre- 
sponds to the well known period T ~ 88y [23]. Compar- 
ing this with the Fig. 1 one can conclude that the period 
(1933y-...) of the solar hyperactivity is close to its end. 
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Fig. 8: Spectrum of the sunspots number fluctuations in the 
In- linear scales (the reconstructed data for the last 11,000 years 
have been taken from [20]). The straight line is drawn to indi- 
cate the exponential law Eq. (3). 
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